#include <iostream>
#include<windows.h>
#include<stdlib.h>
using namespace std;

const int N=5000;
double b[N][N],a[N],sum=0.0;

void init()
{
    for(int i=0;i<N;i++)
    {   a[i]=i;
//        for(int j=0;j<N;j++)
//        {
//            b[i][j]=i+j;
//        }
    }

}

int main()
{
    long long head, tail , freq ; // timers
    init ();
    // similar to CLOCKS_PER_SEC
    QueryPerformanceFrequency((LARGE_INTEGER *)&freq );
    // start time
    QueryPerformanceCounter((LARGE_INTEGER *)&head);
    double sum1=0.0,sum2=0.0,sum3=0.0,sum4=0.0;
for(int n=0;n<10;n++)
{

    sum1=0.0,sum2=0.0,sum3=0.0,sum4=0.0;
    for(int i=0;i<N;i+=4)
    {
        sum1+=a[i];
        sum2+=a[i+1];
        sum1+=a[i+2];
        sum2+=a[i+3];
    }
    sum=sum1+sum2+sum3+sum4;
//    double sum1[N],sum2[N];
//    for(int i=0;i<N;i+=2)
//    {
//        sum1[i]=0.0;
//        sum1[i+1]=0.0;
//        sum2[i]=0.0;
//        sum2[i+1]=0.0;
//    }
//    for(int j=0;j<N;j++)
//    {
//        for(int i=0;i<N;i+=2)
//        {
//            sum1[i]+=b[j][i]*a[j];
//            sum2[i]+=b[j][i+1]*a[j];
//        }
//    }
//    for(int i=0;i<N;i+=2)
//    {
//       sum[i]=sum1[i]+sum2[i];
//       sum[i+1]=sum1[i+1]+sum2[i+1];
//    }
}
    // end time
    QueryPerformanceCounter((LARGE_INTEGER *)&tail );
    cout<<"Time:" << ((tail- head) * 1.0 )/ freq << "s" << endl;


    system("pause");
    return 0;
}


